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The spatial clustering of galaxies in galaxy clusters implies that the background of infrared (IR) 
light in the intracluster medium (ICM) may exceed the universal background. Cosmic rays injected 
within the ICM propagate diffusively and at low enough energies are trapped there for cosmological 
times. The photopion production interactions of cosmic rays with the IR photons are responsible 
for the generation of neutrinos whose detection may shed some light on the origin and propagation 
of high energy cosmic rays in the universe. Here we discuss our calculations of the flux of neutrinos 
from single clusters as well as the contribution of photopion production in clusters of galaxies to the 
diffuse neutrino background. 

PACS numbers: 



I. INTRODUCTION 

Clusters of galaxies are peculiar storage rooms of cos- 
mic rays Q, Q: most non-thermal protons and nuclei 
injected in the ICM remain confined there for cosmo- 
logical times, and can either interact with the gas and 
photon backgrounds or be re-energized by episodic phe- 
nomena, such as mergers with other clusters 0, Q . The 
spatial clustering of 10-1000 galaxies within the cluster 
volume also enhances the amount of infrared and opti- 
cal light in that region, thereby increasing the probabil- 
ity of interactions of charged particles with these pho- 
tons. The enhanced cosmic ray and photon densities 
make clusters of galaxies the ideal laboratory for the gen- 
eration of high energy radiations and neutrinos. The de- 
tection of a neutrino flux from clusters of galaxies would 
be a precious tool to weigh the unknown non-thermal 
content of these large scale structures. As discussed in 
Sec. [n] the amount of cosmic rays trapped in the intr- 
acluster medium is still very uncertain, and only weakly 
constrained by present observations of hard X-rays and 
gamma rays, most likely generated by high energy elec- 
trons 0, [||. The detection of neutrinos on the other 
hand would allow us to have an independent handle to 
infer the energetic non-thermal content which is in the 
form of hadronic cosmic rays. Gamma rays in the 1-1000 
GeV range are also expected to carry precious informa- 
tion about the amount of cosmic rays trapped in the ICM 

Hi 

In this paper we investigate the most promising mecha- 
nism for the generation of high energy neutrinos, namely 
the production and decay of charged pions in proton- 
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photon inelastic interactions. The relevant photon back- 
ground in the ICM is provided by the IR and optical light 
emitted by galaxies within the cluster. The maximum 
flux of neutrinos from inelastic pp interactions allowed 
by the observations of the diffuse gamma ray background 
was calculated in 0, Q . The calculation of the neutrino 
flux is carried out both analytically and through a Monte 
Carlo simulation of the diffusive motion of high energy 
particles in the turbulent magnetic field of a cluster of 
galaxies. Such propagation was also investigated in [t| 
with the purpose of determining the flux of gamma and 
hard X-rays generated as secondaries of the high energy 
protons. 

The paper is organized as follows: in Sec. [n]we review 
the current knowledge of the cosmic ray propagation and 
energetics in the ICM. In Sec. IHII we illustrate our calcu- 
lations of the IR light background and its spatial distri- 
bution in the volume of the cluster. In Sec. HVI we outline 
a benchmark analytical calculation of the neutrino flux, 
in order to emphasize the dependence of the result upon 
poorly known quantities, such as the diffusion coefficient 
and the CR energy content of the cluster. A Monte Carlo 
calculation of these neutrino fluxes is described in Sec. E] 
In Sec. I VII we calculate the maximum expected contribu- 
tion of these processes of neutrino production in clusters 
of galaxies to the diffuse extragalactic neutrino flux. We 
conclude in Sec. IVIII 



II. COSMIC RAYS IN THE INTRACLUSTER 
MEDIUM 

The presence of a non-thermal component in clusters 
of galaxies is best shown by the observation of regions 
of extended radio emission, called radio halos, with size 
1-2 Mpc, often wider than the X-ray emitting region (see 
|10| for a recent review). The spectrum of the radio emis- 
sion is typically a power law with a high frequency cutoff 
in the 1-10 GHz range and it results from synchrotron en- 
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ergy losses of high energy electrons in the magnetic field 
of the ICM. Measurements of the rotation measure pro- 
vide magnetic fields which are of the order of several /zG 
[Til IT^ | , although several of these inferred values depend 
on assumptions made on the structure of the magnetic 
field and on the density of background electrons. In a few 
clusters, the radio emission is accompanied by X-rays in 
excess of those of thermal origin (due to bremmstrahhmg 
emission). These non-thermal X-rays are most likely the 
result of inverse Compton scattering (ICS) of the same 
electron population responsible for the radio emission. 
The main photon target is represented by the photons 
of the cosmic microwave background (CMB), while the 
infrared light, even in the case of the enhanced fluxes 
considered here, does not play a significant role. In the 
few cases in which both a radio halo emission and a hard 
X-ray emission are detected, an estimate of the magnetic 
field can be obtained ^| • Typical values are in the range 
0.1-1 fiG. These fields should be considered as averages 
over the volume of the cluster. Higher fields are likely to 
be present in the central regions of clusters. 

While non-thermal electrons are easily visible in clus- 
ters, a much harder task is to identify hadronic cosmic 
rays. The smoking gun showing their presence, namely 
a gamma ray emission associated with the decay of neu- 
tral pions, is still lacking. This absence of detection is 
however still compatible with most models of injection 
of cosmic ray protons in the ICM: the expected gamma 
ray emission becomes comparable with the EGRET sen- 
sitivity only if the protons in clusters have a total energy 
comparable with the thermal energy in the virialized gas. 
Lower fluxes should however be detectable with the up- 
coming GLAST satellite (0 E3, see also @ for a recent 
review). One of the main reasons for the interest in the 
cosmic ray population in clusters is related to the fact 
that the presence of magnetic fields is able to confine cos- 
mic rays for times in excess of the age of the cluster: these 
large scale structures work then as storage rooms for cos- 
mic rays 0, 0] • This implies that the energy density of 
cosmic rays in clusters increases with cosmic time. The 
potential sources of cosmic rays in clusters have been con- 
sidered in [| and the corresponding energy densities were 
estimated. The efficiency of the confinement was also in- 
vestigated: the maximum energy for which the confine- 
ment is effective, namely the confinement time exceeds 
the age of the cluster (roughly the age of the universe), 
is strongly dependent upon the energy dependence of the 
diffusion coefficient. For Bohm diffusion and fiG fields, 
the confinement remains effective up to extremely high 
energies, while a Kolmogorov spectrum of magnetic fluc- 
tuations leads to confinement only up to energies of the 
order of several TeV. 

Normal galaxies in clusters are expected to contribute 
a total energy in the form of confined cosmic rays e ~ 
LpA^gaiio, where A^ai ~ 100 is the typical number of 
galaxies in a cluster, L p ~ 3 x 10 40 erg s _1 is the cos- 
mic ray luminosity of a galaxy like ours. We then obtain 
e w 10 60 erg, more than three orders of magnitude less 



than the total energy in the thermal gas. A single active 
galaxy in a cluster should easily contribute a cosmic ray 
luminosity of ~ 10 44 erg s~ x . At any time it is reasonable 
to assume that at least one of the galaxies in the cluster 
is in an active phase, so that the total energy contributed 
by active galaxies is expected to be ~ 3 x 10 61 erg. A sim- 
ilar estimate is obtained if cosmic rays are accelerated 
through first order Fermi acceleration at the accretion 
shock in the outskirts of the cluster, if the acceleration 
efficiency is ~ 10%. These shocks form in the outskirts 
of clusters and are the result of the propagation of the 
information on the virialization of the core towards the 
outer regions . Accretion shocks propagate in the cold 
intergalactic gas and therefore their Mach number can be 
very high. In the context of the linear theory of parti- 
cle acceleration at non-relativistic shocks, the expected 
spectrum of accelerated particles in the limit of very large 
Mach numbers is a power law with slope a « 2. The cos- 
mic ray induced modification of the shock, particularly 
important for high Mach numbers, may gene rate higher 
efficiencies and non-power- law spectra 0, H3, UK • 

Shock waves are also expected to appear due to the su- 
personic infall of two or more clusters in a merger event. 
The typical strength of these shocks is much lower than 
for accretion shocks, and the expected spectra of accel- 
erated particles is steep 0, |2(]]. The energy liberated 
in the gravitational form during these events is of order 
~ GMg u /R vlI ~ 10 64 erg, which is in fact of the same 
order of magnitude of the thermal energy in the cluster, 
as should be expected since these mergers are believed 
to be responsible for the heating of the ICM to its virial 
value. Some fraction of this energy may be in the form 
of accelerated protons, but the efficiency for particle ac- 
celeration is expected to be rather small because of the 
low Mach numbers, which are also responsible for steep 
spectra 

From the observational point of view, the most severe 
constraint to the energy density of cosmic rays, though 
model dependent, is represented by the radio emissivity 
at a few GHz frequency. Radio emission at this frequency 
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is generated by electrons with energy E e — 37GeVB M , 
where is the magnetic field in /iG. Observations of 
the Coma cluster show a cutoff in the volume integrated 
radio spectrum, at v ~ 5 GHz [2l| . If cosmic rays are 
present, a contribution to the radiating electrons comes 
from secondary electrons, generated from pp inelastic col- 
lisions. The spectrum of secondary electrons however has 
no cutoff at such low energies, therefore the radio emis- 
sivity at high frequencies imposes a rather severe con- 
straint on the amount of secondary electrons, and there- 
fore of primary cosmic rays in the region 10 — 1000 GeV 
[22| . This bound implies that cosmic rays amount to less 
than (1 — 10)% of the thermal energy in the cluster, de- 
pending upon assumptions on the spectrum at injection 
and the morphology of the magnetic field in the cluster. 
These limits are slightly tighter if reacceleration of the 
secondary electron-positron pairs is taken into account 
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III. THE INFRARED LIGHT IN THE 
INTRACLUSTER MEDIUM 

In this section we illustrate our determination of the 
InfraRed (IR) background in the intracluster medium. 
We start with adopting the Spectral Energy Distribu- 
tion (SED) of individual galaxies following [24]]. We then 
carry out the convolution of these SEDs with the spatial 
distribution of galaxies in a sample cluster of galaxies. 

The SED of a 10 n £o normal galaxy as obtained by 
[Hill is illustrated in Fig.U The authors of Ref. 0HI 
state that the SED they find is valid only at wavelengths 
longer than 4/Ltm [2 II . but after comparing it with other 
models and data we decided to use it at longer wave- 
lengths as well, since the differences were found to be 
smaller than any other uncertainty of astrophysical ori- 
gin. The same comment holds for the differences due to 
different types of galaxies in their model. 

We used lfim as upper energy for the background ra- 
diation since at higher energies its spectrum is much 
steeper and even with a steep proton spectrum there is 
no appreciable contribution to the neutrino production. 
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FIG. 1: Model spectrum of a normal galaxy of luminosity 
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In order to determine the IR background within the 
cluster of galaxies, we calculate the convolution of the 
given SED from our prototypical galaxy with the spatial 
distribution of galaxies. Following 9] we assume that the 
galaxy distribution follows closely the distribution of gas, 
for which we use: 

r/ r \0.51/ r \ 0.72 / rN 0.58-1-1 

with n = lOkpc, r 2 = 250 kpc and r 3 = IMpc. Fig. H 
shows the energy density at the peak of the SED as a 
function of the distance from the center of the galaxy 
cluster, in the assumption that N g ~ 1000 galaxies are 
in the cluster. The function shown in Fig. [21 can be used 
to construct a useful model of the IR density as a function 
of the distance to the center of the cluster by integration 
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FIG. 2: Energy density at the FIR peak for a particular re- 
alization (points) and its fit (line) that can be used for other 
realizations. 

over different distance ranges. This is shown in Fig. 
where the sum of the IR emission from cluster sources is 
added to the extragalactic IR background and plotted in 
six circular annuli inside the cluster. 
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FIG. 3: Energy density of the IR (cluster + extragalactic) 
for six zones inside the cluster. Points show the extragalactic 
background from Ref. |2rj| . 



IV. NEUTRINO PRODUCTION: AN 
ANALYTICAL CALCULATION 

In this section we present an analytical estimate of the 
neutrino flux generated as a result of the photopion in- 
teractions of high energy protons with the IR and CMB 
photons in the ICM. For simplicity we assume that there 
is only one dominant source of energetic protons in the 
center of the cluster and that the IR photon background 
is constant within the inner 500 kpc of the cluster. We 
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also assume that the diffusion coefficient within this re- 
gion is spatially constant. 

Under these assumptions the neutrino production in 
the cluster can be calculated as: 



AE„ 



xn p (E p ,r)——4Trr\ (2) 



where An v jAE v is the average spectrum of neutrinos pro- 
duced by a single proton interaction on the photon back- 
ground, n p (E p , r) is the number of protons at distance r 
from the center and c/£(E p ) is the interaction probability 
per unit time. 

The number density of protons at distance r from the 
central source, n p (E,r), can be easily determined by 
solving the diffusion equation 0, obtaining: 



Q(E) 
ATrrD(E) ' 



(3) 



This solution is valid for particles that in the time to 
manage to diffuse out to a distance r from the center 
without suffering appreciable energy losses. This means 
particles whose energies satisfy r 2 /(AD(E)) < to- At 
distances larger than \J 'AD(E)to the density of particles 
with energy E drops exponentially. 

Here we assume that the diffusion coefficient has the 
form found in [27j : 



D{E) = 



1/3 



(4) 



where £>* = jr-^(E^)c, r^(E^) is the Larmor radius of 
particles with energy E* and the reference energy E* is 
found by requiring that r^(E^) = L c /5, L c being the 
coherence scale of the magnetic field. The three terms 
in Eq. 01 correspond to different regimes of propagation 
(Kolmogorov, Bohm and pitch angle scattering) depend- 
ing upon the comparison of the Larmor radius of the 
particle and the coherence scale of the magnetic field. 

The upper limit in the integral over the ra- 
dial coordinate in Eq. El is defined as r<iig(E) — 
min(y/ 4D(E)t ,r max ), with r max = 500 kpc and t the 
age of the cluster, comparable with the age of the uni- 
verse, t ~ 10 10 yr. It is easy to show, following that, 
for cosmic ray energies for which confinement works, the 
spectrum of secondaries (gamma rays and neutrinos) is 
not affected by diffusion and is therefore not steepened by 
the energy dependence of D(E). For unconfined cosmic 
rays, diffusion plays a role and the spectrum of secon- 
daries suffers a steepening. 

For a given diffusion coefficient, the minimum energy 
of the particles for which the escape time does not exceed 
to can be estimated by requiring that r^ ax /4Z?(_E) < to- 

2 

In Fig.0]we plot the diffusion time, t(E) = , out of 

a sphere with radius r max as a function of energy, using as 
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FIG. 4: Diffusion time as a function of the proton energy. 
In the energy region considered here, the gray region in the 
plot, the diffusion time is always shorter than the age of the 
cluster. 



a diffusion coefficient the one reported in Eq. The two 
curves correspond to magnetic fields in the cluster Bo = 
1 /iG (lower curve) and B = 5 /iG (higher curve). As it 
is clear from the plot, due to the large energy threshold 
for photopion reaction, all the particles of interest for 
our problem are not confined on time scales of the order 
of the age of the cluster. It is however worth keeping in 
mind that this estimate is rather simplified and that more 
realistic distributions of magnetic field may reflect into 
longer confinement times as we actually find in Sec. 

We used the SOPHIA event generator [2^ to calculate 
the average spectrum of neutrinos produced in a single 
proton interaction on the cluster photon background and 
we parametrized the results as: 



&E U 



( Ep ,E v ) 



45 
E 



■ exp 



\E p /20) 



1.2- 



(5) 



for E v < E p and otherwise. In order to calculate the 
interaction probability per unit time we use the inter- 
action lengths as calculated in the next section and we 
approximate the interaction length in the region with 



3 • 10 14 (^) 12 Mpc, with 



r < 500 kpc as: £(E P ) 
E = lGeV (see Fig.EJ. 

The integration limits for the proton energy in Eq. 
El are: E™ ax = 10 11 GeV and E™ in = max(E^E th = 
10 8 GeV) (for E p < 10 s GeV the neutrino production is 
negligible) . 

For the purpose of numerical calculations, we assume, 
following 0, |2JJ , that the cosmic ray source has a lumi- 
nosity, above energy E = 1 GeV, L p = 3 x 10 44 ergs _1 
(see also Sec. ILTll . 

The neutrino spectrum from Eq. [21 for a = 2.1 is plot- 
ted in Fig. [3 for magnetic field 1 ^G (lower curve) and 
5 [iG (upper curve). Assuming that the cluster is located 
at a typical distance of 100 Mpc, our estimate of the neu- 
trino flux at the Earth is of 0.29 km~ 2 yr _1 for magnetic 
field 1 /iG and 0.7 km -2 yr _1 for magnetic field 5 fiG. 

This result varies somewhat with different choices of 
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FIG. 5: Analytical estimate of the neutrino emissivity from 
proton propagation in a cluster, with slope of the injection 
spectrum of protons a = 2.1. Upper curve for Bo = 5/xG, 
lower curve for Bo = 1 /xG. 

the diffusion coefficient, but not in a substantial way. 
Steeper injection spectra of cosmic rays imply smaller 
neutrino fluxes. 



V. NEUTRINO PRODUCTION: A MONTE 
CARLO APPROACH 

The result obtained in the previous section may some- 
what be affected by the adoption of a more detailed distri- 
bution of the IR photons in the ICM, as well as by details 
of the diffusion in a space dependent magnetic field. In 
particular, both the IR light and the magnetic field are 
expected to be higher in the central region of the clus- 
ter. In order to investigate the effect of these factors on 
the predicted neutrino flux, we performed a Monte Carlo 
calculation of the propagation and the related neutrino 
production. 

The photon background is taken as the superposition 
of the IR background, as calculated in Sec. II I II and the 
usual CMB. We divide the cluster in six concentric zones 
with radii: 100 kpc, 300 kpc, 500 kpc, IMpc, 2Mpc and 
3 Mpc. In each zone we assume the infrared background 
to be constant and equal to the average one in that zone. 
The interaction lengths calculated in each zone are plot- 
ted in Fig. HJ1 

The simulation of the propagation of charged particles 
in the ICM requires a specific choice of the strength, spa- 
tial profile and disordered component of the magnetic 
field. Inspired by the results of [23, El 011 the cluster 
Abell 119, we assume that the total magnetic field scales 
with distance from the center of the cluster as given by 
the flux freezing condition with an electron thermal com- 
ponent modelled as a /3-model: 

B(r)< x (i + - S ) (6) 

with r c = 378 kpc. This result agrees quite well with 
both the theoretical expectation of B{f) cx n (r) 2 / 3 for a 



magnetic field frozen in matter (n e (r) is the electron den- 
sity) and with the simulation results of Ref. (see their 
Fig. 1). As a result, the average magnetic field strength 
in the six zones of the cluster changes from about 4.5 /xG 
in the inner zone to about 0.18 fxG in the outer zone in 
the case of 5 /xG strength in the center. The average field 
values are scaled down by a factor 5 for the 1 /xG case. 
The simulation code does not use these averages, rather 
the magnetic field at the particle location. 

We implemented a turbulent magnetic field with the 
above radial dependence dividing the ICM in cubes of 
50 kpc side, each filled with a random magnetic field of 
average strength (B) = B satisfying a Kolmogorov spec- 
trum with three logarithmic sub-scales (12.5, 25, 50 kpc). 
Given the proton position we calculate the turbulent 
magnetic field in the cube containing the particle and 
then we scale its magnetic field according to Eq. © . For 
more informations and details about the actual imple- 
mentation see Appendix B of Ref. [33 ■ 

We calculate the neutrino yield, Y(E p , E v ), namely the 
spectrum of neutrinos generated during the propagation 
of protons of energy E p in the ICM. This calculation is 
carried out by injecting protons in 30 logarithmic energy 
bins between 10 17 eV and 10 20 eV. We inject 10,000 pro- 
tons for each energy bin. We follow their propagation 
in the ICM until their distance from the center exceeds 
3 Mpc or their energy falls below 10 17 eV. Clearly the 
neutrino yield as defined here depends on the specific 
model for injection of the particles. We limit ourselves 
to considering only the case of a source that generates 
cosmic rays in the center of the cluster. 

The propagation is simulated using the Monte-Carlo 
code described in Ref. [12 • The nucleon-photon 
interactions are simulated with the event generator 
SOPHIA 28]. At each interaction we record the sec- 
ondary products and then at the end of the simulation 
we have the spectrum of the neutrinos produced inside 
the cluster by a proton injected with a given energy, 
Y(E p ,E v ). 

We ran our simulation for two values of the mag- 
netic field in the center of the cluster, Bo = 1 /xG and 
Bo = 5 /xG. In Fig. [fj] we plot the pathlength travelled 
by particles propagating in the cluster as a function of 
energy for the case of injection in the center and for 
Bo = 5 /xG. For each energy bin we calculate the av- 
erage pathlength travelled. In this average we include 
particles that exit the cluster and particles that go below 
threshold. This considerably reduces the pathlengths at 
low energy (below 10 18 eV). 

Since at about ~ 10 17 eV the interaction time becomes 
larger than the age of the universe (see Fig. |f)J) we de- 
cided to use only the protons with energy above 10 17 eV 
to calculate the neutrino flux. Lower energy protons, 
even with a steep spectrum, do not change the predicted 
neutrino flux because at these energies they lose their 
energy mainly due to adiabatic energy losses. 

Fig. [7| shows the fraction of the proton energy that is 
converted to neutrinos before leaving the cluster. For en- 
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FIG. 6: Interaction lengths of protons propagating inside the six concentric zones in the cluster. The solid lines correspond, 
from bottom to top, to the zones with radii: 100 kpc, 300 kpc, 500 kpc, 1 Mpc, 2 Mpc, 3 Mpc. The dashed line is the interaction 
length on the CMB. The dotted line represents the size of the universe. The dot-dashed line represents the approximate 
interaction length we used in the analytical calculation. Data points show the containment distance in the cluster for central 
magnetic field of 5/iG. 
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FIG. 7: The fraction of the proton energy converted to neu- 
trinos inside the cluster is shown with solid lines for 5/j.G 
and 1/j.G central fields. The dotted curve shows the energy 
in neutrinos for propagation in the CMB on 200 Mpc. The 
dash-dotted and dashed lines show the fraction of energy in 
neutrinos for propagation in the cosmic IR background on 
scales of 200 Mpc and 1 Gpc. 



ergies around 2x 10 17 eV this fraction is 0.1% and almost 
linearly increases with \nE p to about 1% at 10 18 eV and 
stays roughly constant at higher energy. The proton en- 
ergy loss in the cluster is a function of the hadronic cross 
section and of the containment time in the cluster. Be- 
cause protons of energy approaching 10 20 eV propagate 
almost in straight lines inside the cluster they are not 
expected to be very efficient in generating neutrinos. 

For the case of protons injected at the center of the 



cluster, the neutrino spectrum is given by 

ji{E v ) = jY l (E pi E„)Q p (E p )dE p , (7) 

where: Y l (E pi E l/ ) is the neutrino yield, for neutrinos of 
type i, at given proton energy and Q p (E p ) is the injection 
rate of protons in the cluster. 

The results of our Monte Carlo calculations are shown 
in Fig. |SJ The upper (lower) panel refers to the case 
Bq = 1 fiG (Bo = 5/iG). Both panels show the well 
known pronounced two peaks |3^| due to the direct neu- 
trino production (higher energy peak) and neutrino pro- 
duction following neutron decay (lower energy peak). 
The v e from neutrons decaying outside the cluster are 
not included. The fluxes of neutrinos of different flavors 
are superimposed and the thick histogram represents the 
total flux. The continuous line is the result of the ana- 
lytical calculations discussed in Sec. IIVI The agreement 
between the analytical and Monte Carlo calculations is 
remarkable, although the former did not include the de- 
cay of neutrons, and therefore has no low energy peak. 

The integrated flux of neutrinos with energy E > 10 5 
GeV as calculated from the Monte Carlo for a cluster 
at a typical distance of 100 Mpc is 0.9 km _2 yr _1 for 
B = 5 [iG and 0.2 km~ 2 yr _1 for B = 1 [iG. One should 
keep in mind that this flux is obtained for a source lu- 
minosity in the form of protons L p = 3 • 10 44 ergs _1 . 
The neutrino flux scales linearly with L p , and may be 
appreciably higher than those predicted here if a sub- 
stantially more luminous source happens to be located 
in a nearby cluster. On the other hand, the neutrino flux 
is dominated by neutrinos with energies around 10 17 eV, 
generated by particles that travel within the ICM for a 
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FIG. 8: Flux of neutrinos from a cluster located at a typical 
distance of 100 Mpc from the Earth. The upper (lower) panel 
refers to the case Bo = 1/iG (Bo = 5/iG). The different 
histograms show the flux of neutrinos of different flavors, while 
the thick histogram shows the total flux. The continuous line 
is the result of the analytical estimate in Sec. IIVI 



substantial fraction of the age of the cluster. The lumi- 
nosity reported above should therefore be interpreted as 
the proton luminosity averaged over a period comparable 
with the age of the cluster and it is unlikely to find values 
of L p much higher than that used here. 



VI. THE DIFFUSE NEUTRINO FLUX 

In this section we illustrate our calculations of the con- 
tribution of photopion production in clusters of galaxies 
to the diffuse neutrino flux. For this purpose, we use the 
results of our Monte Carlo calculations (Sec. [Vj to ob- 
tain a template spectrum to be convolved with the mass 
distribution of clusters of galaxies. For the integral mass 
function of clusters of galaxies we use the following ex- 
pression |3^ |: 



n(> M) —Ax 10" b — 



M 
M 



exp 



-wJ h Mpc 



where M is the total (gravitational) mass of a cluster 
within 1.5 h- 1 Mpc and M* = (1.8 ± 0.3) x 10 14 h^M® 
is a reference cluster mass. This distribution covers the 
range from M = lO 12 ^ 1 !/© to M = 5 x 1O 15 /i _1 M . 
For our numerical calculations we adopt h — 0.7 for the 
dimensionless Hubble constant. The flux of diffuse neu- 
trino radiation from clusters of galaxies, as due to pho- 
topion production can be written as 



W )^/dM^ W ,M), 



(9) 



where we neglect the effect of rcdshift since the uncer- 
tainties in the determination of the neutrino emissivity, 
J„, are very large compared with the effects of cosmology. 
Rather than presenting several cases of diffuse neutrino 
fluxes, depending on the choice of the several parame- 
ters that enter the calculations, we decided here to show 
the maximum neutrino flux. In order to estimate this 
flux, we assume that the maximum allowed cosmic ray 
luminosity in the cluster is 



GM 2 
R v t c \ 



= ••-! - 10" I St 1 



0.1 



M 



5/3 



ergs 



where we used the definition of virial radius as 



3M 



47rA c Sl M /Oc 



1/3 



(10) 



(11) 



with A c w 200, £!m = 0.3 for the matter fraction in the 
universe and p cr = 1.88 x 10~ 29 h 2 gem -3 for the critical 
density. One should keep in mind that luminosities of this 
order of magnitude can only be reached in the context of 
major mergers of clusters of galaxies, namely when two 
clusters with comparable masses merge to form a new 
more massive cluster. Aside from the obvious fact that 
these events are rather rare, one should also remember 
that the shock waves that develop during these events 
and that are responsible for both the heating of the intr- 
acluster gas and the possible acceleration of cosmic rays, 
have relatively low Mach numbers 0] (e.g. M ~ y/2— 2). 
The spectrum of particles accelerated at these shocks is 
typically very steep , and certainly quite steeper that 
E~ 2A , used in our calculations of the neutrino emissiv- 
ity. Eq. E3is therefore expected to be an absolute up- 
per limit to the cosmic ray luminosity (or equivalently to 
the cosmic ray energy density averaged over the age of 
the cluster, t c \ w 10 10 yr), in particular if used together 
with relatively flat spectra (e.g. oc E~ 2A ). These large 
energy densities would also be in contradiction with the 
upper bounds to cosmic ray pressure found in 22] for the 
Coma cluster, although these limits are somewhat model 
dependent. Lower cosmic ray luminosities and flatter in- 
jection spectra are likely to be obtained if particle accel- 
eration takes place at shock waves that develop at the 
outskirts of clusters 0, 0, |3j| , due to their larger Mach 
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FIG. 9: Upper limit to the neutrino flux from photopion pro- 
duction in clusters of galaxies for maximum energy of the 
protons E max = 10 20 eV (thick histogram) and E max = 10 19 
eV (thin histogram). We also show the flux of atmospheric 
neutrinos (hatched region) including the contribution from 
prompt neutrinos (from |3g). The sensitivities of AMANDA 
II and ICECUBE are also shown (from [37|). 



numbers (namely lower temperatures of the background 
unshocked gas). 

The estimate of our upper limit on the diffuse neutrino 
flux from photopion production in clusters of galaxies is 
shown in Fig. [5] the thick (thin) histogram represents 
our upper limit when the maximum energy of protons in 
clusters is E max = 10 20 eV (£ max = 10 19 cV). We also 
show there the upper bound on the neutrino flux from pp 
interactions as found in fil (dotted line). Note however 
that the upper limit in [lj was found by requiring that 
gamma rays generated in proton-proton inelastic scat- 
terings in the ICM would saturate the EGRET diffuse 
gamma ray background. This would require £ 0.3 in 
Eq. I1UI and injection spectrum E~ 2A . As we already 
stressed, this situation is totally unphysical, as also dis- 
cussed in but represents a solid upper limit to the 
flux of neutrinos from clusters. 

In Fig. Ejth e flux of atmospheric neutrinos is also plot- 
ted (from 36] ) in order to emphasize the flatness of the 
neutrino spectrum as predicted here and the energy re- 
gion in which neutrinos from photopion production in 
clusters may overcome the atmospheric neutrino flux (the 
dark-hatched region is the standard atmospheric neu- 
trino flux, while the light-hatched region refers to prompt 
neutrino emission from charmed mesons). This happens 
at energy E > 5 x 10 6 GeV. The upper limits that we 
show are still above the ICECUBE sensitivity (but below 
AMANDA II), as shown in the plot. 

In obtaining the upper limit neutrino flux in Fig. |5J we 
did not consider however another important constraint: 
as we stressed in Sec. IIVI there are cosmic rays that are 
able to escape from the cluster volume in times shorter 
than the age of the cluster. Following Fig. El we adopt 
as a fiducial energy for the escape 10 18 eV: particles with 




E[eV] 



FIG. 10: Flux of cosmic rays escaping clusters for E mElx = 
10 20 eV and E 

max — 10 19 eV. The curves are normalized to 
the data of AGASA and Akeno. The factors 40 and 14 are 
the suppression factors needed to avoid overproduction of the 
extragalactic cosmic rays. 



larger energies are not trapped inside the cluster and con- 
tribute to the flux of extragalactic cosmic rays at the 
Earth. It is important to realize that the spectrum of 
the escaping particles is very close to the injection spec- 
trum at the source and does not resemble the spectrum 
of cosmic rays in the ICM (Eq. |3J • The effect of energy 
losses on the spectrum of escaping cosmic rays is neg- 
ligible. We calculated the flux of cosmic rays escaping 
clusters and reaching the Earth following the well known 
approach to propagation illustrated in |38| , which allows 
us to take into account photopion production, Bethe- 
Heitler pair production and adiabatic losses on cosmo- 
logical scales. The results of our calculations are plotted 
in Fig. 1 101 for a maximum energy of the primary protons 
10 20 eV and 10 19 eV. The low energy (~ 10 18 eV) cutoff 
in the spectra reflects the lack of confinement of higher 
energy cosmic rays. The curves are normalized to the 
data of AGASA and Akeno: in order to avoid exceeding 
the observed fluxes of cosmic rays at the Earth, the flux 
from clusters needs to be suppressed by a factor ~ 40 for 
£ max = 10 20 and by a factor - 14 for E max = 10 19 . The 
neutrino fluxes in Fig. [5] need to be lowered by the same 
factors, so that the predicted neutrino fluxes from pho- 
topion production in clusters drop below the sensitivity 
limit for ICECUBE. 



VII. DISCUSSION AND CONCLUSIONS 

We calculated the neutrino flux that could be expected 
from the photopion interactions of very high energy pro- 
tons with the infrared light and the cosmic microwave 
background in clusters of galaxies. The main target for 
these interactions is represented by the infrared photons 
generated by galaxies in the cluster itself. This back- 
ground can be much higher than the cosmic infrared 
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background. This argument, added to the effective con- 
finement of cosmic rays within the ICM for long times 
make clusters interesting targets for neutrino telescopes. 

The calculations have been carried out both analyti- 
cally and with a Monte Carlo code for the propagation 
and interaction of high energy protons. There are many 
uncertainties involved in these calculations: the cosmic 
ray luminosity in the ICM averaged over cosmological 
times is very uncertain, though somewhat limited by mul- 
tifrequency observations of some clusters. These obser- 
vations suggest that less than 1 — 10% of the thermal 
energy content of cluster is in the form of cosmic rays. 
This limit is rather weakly dependent on the diffusion co- 
efficient, which instead strongly affects the escape times 
out of clusters. 

The confinement time is important for our calculations 
because neutrinos are mainly produced by cosmic rays 
with energy in excess of 10 17 eV, where the confinement 
time is expected to be comparable with or shorter than 
the age of a cluster (Fig. 0). The extent to which this 
is true depends on the choice of the diffusion coefficient 
which is unknown and poorly constrained due to our ig- 
norance of the strength, structure and coherence scale of 
the magnetic field in the ICM. As shown in Fig-El for rea- 
sonable assumptions on this magnetic field, the diffusion 
coefficient in clusters is expected to have a dependence 
on momentum which is roughly linear in the energy re- 
gion from 3 x 10 18 to 3 x 10 19 eV. We have checked the 
correctness of this finding in toy propagation models and 
concluded that without particle energy losses the con- 
tainment time is exactly E^ 1 , despite the fact that the 
model of magnetic field turbulence adopted for the scat- 
tering of the particles is of Kolmogorov type. While we 
believe that we describe well the proton propagation in 
the transition region, the diffusive part of the propaga- 
tion might be not fully realistic. If we overestimated the 
diffusion time of cosmic rays out of the clusters, then the 
neutrino flux is correspondingly reduced. This applies 
equally well to all neutrino flavors. 

The location where cosmic rays are injected also 
changes the predicted neutrino fluxes but only in a rather 
marginal way. If the sources are spread over the cluster 
volume rather than being concentrated in the center as 
we assumed, then the neutrino fluxes are estimated to 



be a factor ~ 2 lower than predicted here. On the other 
hand it is reasonable to think that the IR background 
could have been higher in the past epochs due to lumi- 
nosity evolution of the galaxies in clusters. This evolu- 
tion might enhance the rate of neutrino production in 
the past, that would reflect into a higher diffuse neutrino 
flux for a given cosmic ray luminosity. 

The neutrino fluxes from single galaxy clusters as cal- 
culated here were found to be undetectable with the cur- 
rently planned and constructed detectors, such as Ice- 
Cube [3!j and the European km 3 detector |4(| . Their de- 
tection would require the development of detection tech- 
niques able to provide much higher detection volumes, 
such as the radio detection technique |41j . 

The superposition of the tenuous neutrino fluxes from 
all clusters in the universe generates a diffuse neutrino 
flux and turns out to be somewhat more interesting that 
the single sources. We estimated this diffuse neutrino 
flux in Sec. lVII and our results are plotted in Fig. |5J The 
upper limit to the neutrino flux appears to be detectable 
by IceCube. On the other hand, the cosmic rays that are 
responsible for the production of these neutrinos also es- 
cape the cluster and may reach us as cosmic rays. When 
the flux of these cosmic rays is calculated, it exceeds the 
observed fluxes of cosmic rays by a factor ~ 40 if the 
maximum energy of protons is 10 20 eV and by a factor 
~ 14 if this maximum energy is 10 19 eV. The neutrino 
fluxes are expected to be suppressed by the same fac- 
tors, unless some other process inhibits the propagation 
of these cosmic rays on cosmological scales. One example 
of such processes, although not a viable one, could be pro- 
vided by the presence of a non negligible magnetic field 
in the intergalactic medium (spread over scales compara- 
ble with the loss length of particles with energies ~ 10 19 
eV). However, the value of the magnetic field needed to 
shrink the so-called magnetic horizon of particles of en- 
ergy ~ 10 19 eV to, say, 1/10 of the size of the universe is 
of the order of ~ 10~ 8 G, which appears to be unreason- 
ably large for the magnetic field averaged over the size of 
the universe. 
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